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Why the braking indices of young pulsars are less than 3? 
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Abstract. In this letter we discuss two possible reasons which cause the observed braking indices n of young radio 
pulsars to be smaller than 3: (a) the evolving spin-down model of the magnetic field component B± increases with 
time; (b) the extrinsic braking torque model in which the tidal torques exerted on the pulsar by the fallback disk, 
and carries away the spin angular momentum from the pulsar. Based on some simple assumptions, we derive the 
expression of the braking indices, and calculate the spin-down evolutionary tracks of pulsars for different input 
parameters. 

Key words, pulsars: general — Stars: magnetic fields — Stars: rotation 



^ ! lll968l) 



Oh 
i 

o 

CO 



X 



1. Introduction 

Since the first radio pulsar was discovered bv lHewish et alJ 
1968|h continued radio observations at several wave- 
lengths have not only increased the number of detected 
pulsars but also amassed a wealth of details about the 
temporal, spectral and polarization properties of the ob- 
served pulses. It was then quickly established that ra- 
dio pulsars are rapidly rotating strongly magnetized neu- 
tron stars, powered by their rotational kinetic e nergy 
and losing energy by magnetic dipole radiation l|Goldl 
Il968(l . The spin angular velocity of pulsars decreases with 
time, and the spin-down m ay be written by the relation 
l|Manchester fc Tavlodli977|) . 



-Kfl n , 



(1) 



where f2 is the spin angular velocity, K is a positive pa- 
rameter depending on the magnetic dipole moment, and 
n is the braking index. The braking index n can be a mea- 
sured for young radio pulsars by differentiating equation 
(1), 



(2) 



If the spin-down torque entirely results from mag- 
netic dipole radiation with a constant magnetic field, 
the predicted braking index is n = 3. However, if the 
braking is due to a re l ativis tic stellar wind, the braking 
index n = 1 l|Michell Il969l) . The observed braking in- 
dex n, which has been determined for several young ra- 
dio pulsars, are all less than 3. For example, the brak- 
ing in dex of the Crab pulsar PSR 0531+21 is 2.509 ± 
0.0 01 llLvne et al.lll993lh n= 2.140 ± 0.009 for PSR 0540- 
69 dLivingstone et alJl2005h . n=2.837 ± 0.001 for PSR 



1509-58 llKasni et alJll994lh n=2.91 ± 0.05 for PSR 1119- 
6127, n = 1.4 ± 0.2 for the Vela pulsar PSR 0833-45 
ilLvne et alJll99ri . and n = 2.65 ± 0.01 for PSR J1846- 
0258 w hich located at the cent er of supernova remnant 
Kes 75 (|Livingstone et al.ll2006jh 

In previous works, many authors proposed various 
models to explain the discrepancy between observed 
and theoretical predicted value of n. IPeng et alJ l|l982h 
suggested that neutrino and photon radiation com- 
ing from superfluid neutr ons may brake the pulsars. 
iBlandford fc Romanil l)l988(l propose d that n < 3 is due to 
a m ultipole field and f ield evolution. iMenou et al ] il200lft 
and lAloar et alJ l)200l|) presented a model that both mag- 
netic dipole radiation and the propeller torque applied by 
the debris disk formed after supernova ex plosion c an cause 
spin-d own of pulsars. IXu fc Qiacl l|200d) and IWu et all 
(2003) pointed out that the combination of dipole ra- 
diation and the unipo lar generator may cause n < 3. 
lAllen fc Horvathl lll997h suggested that a variation of the 
torque function is important attribution for low braking 
index, and so on. 

In this letter, we suggest that the braking index may 
be less than 3 due to increasing of magnetic field strength 
of pulsars or braking of fallback disks surrounding pulsars. 
In next section, we describe two simple models to explain 
why the observed braking index n < 3. In section 3 we 
present the calculated results for the braking model of the 
fallback disk. In section 4 we make a brief discussion and 
summary. 



2. Model 

Pulsar braking process is still poorly known even more 
than several decades after it was discovered. The best 
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radiation model is the magnetic dipole radiation, one 
which has been proven to be extremely succes sful in 
expla i ning the observed pro perties of pulsars l)Pacini 
Il9fi7t ICnnn fc Ostrikerl ll9fi9). The rot ating pulsars are 
brake d by magnetic dipole radiation I Ostriker fc Guml 
Il969l) or particle outflow fr om the poles of a dipole field 
ijGoldreich fc Jullianl Il969j) which has a constant field 
strength. However, many works propose that there are 
other energy loss mechanisms and braking torques besides 
magnetic dipole radiation, 



m = -Kn 3 + r(fi), 



(3) 



where K = 2B 2 R 6 sin 2 0/3c 3 ; 9 is the inclination of the 
magnetic axis with respect to the rotation axis; /, B, R 
are the momentum of inertia, the surface magnetic field 
strength, and the radius of the pulsar, respectively; c is 
the velocity of light, and T(Q) is other braking torque. 

We suggest that the observed braking index is less than 
3 possibly due to the following two reasons: (a) K — K( fl), 
T(il) = 0, i.e. the other braking torque does not exist 
except for magnetic dipole radiation; (b) K = constant, 

r(n) ^ o. 

2.1. Case A: K = K(Q),T(Q) = 

If the momentum of inertia of the pulsar is a constant, 
by differentiating equation (3), from equation (2) we can 
obtain the observed braking index 



n = 3 



rflnK 
dlnfT 



(4) 



If the second term on the right-hand side of the above 
equation is negative, therefore n < 3 for young pulsars. 

Assuming that the surface dipole magnetic field 
changes with angular velocity due to some unknown phys- 
ical mechanisms, i?, and the braking index n is constant 
during the pulsar spin-down, we can derive the evolution 
of magnetic field component B± = Bsin6 in the direc- 
tion perpendicular to the rotation axis with spin angular 
velocity of pulsar from equation (4) 



cn— 



(5) 



where C is the integrating constant determined by the 
parameters of the pulsar. Setting the initial spin angular 
velocity is Oo, and the magnetic field component is -Bj_,o 
when pulsar was born, then 



B i = B 



n 



(6) 



From equation (6), the magnetic field component B± 
is constant during spin down if n = 3. For young radio pul- 
sars, n < 3, f£ < Slo due to spin-down, we therefore can ob- 
tain B± > B± q. Namely, the magnetic field vertical com- 
ponent B± of young pulsars would increase if the braking 
torque is only magnetic dipole radi ation and the assump- 
tion m entioned in the above is true l|Blandford fc Romanil 
1988). 



2.2. Case B: K = constant, T(Q) / 

Actually, other brakin g torque possib ly exist besides mag- 
netic dipole radiation. Michel (1988) proposed that there 
may be disks sur rounding radio puls ars, produced by su- 
pernova fallback. iMenou et alJ l|200lf) applied the fallback 
disk model and the propeller torque to explain the lower 
braking index of young radio pulsars. Following this idea 
we assume that there is a fallback disk residing outside the 
light cylinder of young pulsars, and tidal torques instead 
traditional propeller torques are put on the fallback disk, 
and carry away angular momentum from the pulsar inside 
it. 

Same as sta ndard thin disk, the viscou s torque in the 
fallback disk is l|Shakura fc Sunvaevl ll973|) 



T 



2itri>Yir— — r, 
dr 



(7) 



where v is the viscosity and £ is the surface density, Q r — 
(GM/r 3 ) 1 / 2 is the local rotational angular velocity in the 
fallback disk, and M is the mass of pulsar. 

If the gravitational interaction of the pulsar with the 
fallback disk only occurs at the inner edge n of the fallback 
disk, angular momentum from pulsar feed into the fallback 
disk at a rate proportional to the surface density £; at the 
inner region of the fallback disk. As a result of continuous 
input of mass at the inner edge of the fallback disk, the 
viscous torque exer ted on the pulsar by the fallback disk 
can be shown to be l)Taam fc Spruitl2(~)0lHSDruit fc Tarim] 
I2QQ3) 



R 



1/2 



CIR M 



1/2 



(8) 



where, t denotes the spin-down age of the pulsar, M is the 
mass inflow rate (assumed to be constant here), i v j is the 
viscous timescale at the inner edge of disk. 

Assuming that the fallback disk is hydrostatically sup- 
ported and geometrically thin with a pressure scale height 
to radi us ratio of/3 = H-Jr\, we estimate the viscos- 
ity v x l|Chen et alJ l2006|) at the inner edge of the fall- 
back disk using standard a prescription presented by 
IShakura fc SunvaevI l)l973|) 



v x = a ss (3 2 nR 



•3/2 r l/2 



(9) 



where ass is the viscosity parameter. In the following cal- 
culations we set ass = 0.001. 

By simple algebra, from the above equations we obtain 
the viscous timescale at the inner edge in the fallback disk, 

■ (10) 



The total braking torque consisting of magnetic dipole 
radiation and the fallback disk braking can be written as 



m = -Kn 3 + r d (fi), 

with 

T d (fi) = - 7 - 1 c" 1 / 4 i? 9 / 4 /3(3as S t/4) 1/2 Mr! 7/4 , 



(11) 
(12) 
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Fig. 1. The pulsar evolution in the P - P diagram for 
different input parameters /3 and the constant mass inflow 
rate M in the fallback disk. 
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Fig. 2. Evolutionary tracks of braking index for different 
input parameters (3 and the constant mass inflow rate M 
in the fallback disk. 



where 7 = (n/Pic) 1 / 4 , R\ c = c/fl is the light-cylinder 
radius of the pulsar. We assume the inner edge in fall- 
back disk to be located at the light-cylinder radius 
l|Menou et al.ll200lj) . i.e. 7=1. 

Similar as in the above subsection, differentiating 
equation (11), from equation (2) we can get 



5 T 

n = 3- (- + 



(13) 



where r = — r2/(2S7) is the canonical spin-down age of 
pulsars, 77 = T d (n)/(~Kn 3 + T d (Q)) is the efficiency of 
the disk braking torque. The second term on the right- 
hand side of the above equation is positive, resulting in 
the braking index less than 3. 



3. Calculated Results for Case B 

We have adopted the semi-analytical method to calcu- 
late the evolution of P, P, and n for radio pulsars with 
a fallback disk. We set some typical values for pulsars, 
B± = i0 12 G, M = I.4.M0, R = 10 6 cm, i.e. K ~ 
2.5 x 10 28 g-cm 2 s I « 10 45 g • cm 2 ; and (3 = 0.02, 0.03, 0.04 
l(Belle et al.l 12004) in the calculations. Setting the initial 
spin period of the pulsar P = 0.02s, we obtain the pulsar 
evolution in the P— P diagram and the evolutionary tracks 
of the braking index in Figs. 1 and 2 for different values 
of M and f3, respectively. We can find that the braking 
index is all the way less than 3, and first increases then 
decreases during the spin-down from Fig. 2. The larger M 
and (3, the larger P and the smaller n. 

For the Crab pulsar, we can derive the current braking 
index since its age is a certain value in the case of a steady 
inflow in fallback disk . Set t — 930yr, the mass inflow rate 
M = 3x 10 17 g-s-\ and = 0.03, 0.04, 0.05, we derive the 
braking index are n — 2.60, 2.49, 2.39 from equation (13), 
which can account for the measured value within 5%. 
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Fig. 3. Evolutionary tracks for four young radio pulsars, 
the solid, dot, dashed-dot, and dashed-dot dot lines corre- 
spond to the Crab pulsar, the Vela pulsar, PSR 0540-69, 
and PSR 1509-58, respectively. The four open stars denote 
the present location of four young pulsars. 



4. Discussion and Summary 

If the radiation process of pulsars is pure magnetic dipole 
radiation, and /, M, R are constant, the change of mag- 
netic field strength is the possible reason of lower brak- 
ing index. From f2 oc f2 n , we can derive the evolutionary 
tracks of four young pulsars in P — P diagram (see Fig. 
3). Anomalous X-ray pulsars (AXPs) and soft Gamma- 
ray repeaters (SGRs) are a small class of pulsars with the 
relatively n arrow period range of 5 — 12s and a higher spin- 
down rates ^Manchester 120011) . The spin-down luminosity 
E = 4tt 2 PP/P 3 of AXPs and SGRs is much less than their 
X-ray luminosity. iThompson fc Duncan! l|l996|) suggested 
that they are powered by decay of strong magnetic fields. 
It is clear that the magnetic field strength of Vela pulsar 
can exceed 10 14 G when P ~ 5 — 6s from Fig. 3, suggesting 
that magnetars as AXPs and SGRs possibly evol ve from 
normal radio pulsars as Vela JLin fc Zhang Il2004h . 
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However, the neutron star may lose spin angular mo- 
mentum not only by magnetic dipole radiation but also by 
other braking mechanisms. Many authors proposed that 
young radio pulsars could be surroun ded by a remnant 
disk genera ted by supernova fallback l)Michel fc Dessleil 
I1981L Il983h . which may contribute the braking torque 
causing spin-down of neutron stars. Assuming that a con- 
stant rate of material feeding into fallback disk, we derive 
the tidal torque exerted on the pulsar by the fallback disk. 
Our calculations show that the braking torque is signifi- 
cantly influenced by the input parameters M and /3. The 
larger values of M or (3 can yield a larger braking torque 
and a smaller braking index n. Fig. 1 shows that large M 
or (3 can cause the generation of pulsars with high spin- 
down rate, like SGRs and AXPs. The evolution of the 
braking index shows that they are all less than 3 during 
spin-down (see Fig. 2). 

If some of the so-called magnetar, i.e. SGRs and AXPs, 
evolve from normal young pulsars as described by case a, 
their surface magnetic field should be strong as iS 10 14 
G, while in case b they may have usual magnetic field 
~ 10 12 G but with high spin-down rates. Measurement of 
cyclotron lines in magnetars may present constrains on 
the two kinds of models. In addition, the presence of fall- 
back disks surrounding pulsars could be detected, since 
a relatively dense fallback disk can produce optical and 
infrared radiation. Optical emission is a common char- 
acter of young pulsars. The Crab pulsar, the Vela pul- 
sar, and PSR 0540-69, which possess low braking index, 
have been shown to pos sess a pulsating optical counter- 
part l|Menou et alJl200llh implying the possible existence 
of the fallback disks. Searching for optical and infrared 
emission from the fallback disks surrounding young pul- 
sars will be of great importance for their evolution. 
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